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Abstract 


The measurement of primary n ± , K 1 , p and p production at mid-rapidity (|y < 0.5) in proton- 
proton collisions at \fs = 7 TeV performed with ALICE (A Large Ion Collider Experiment) at the 
Large Hadron Collider (LHC) is reported. Particle identification is performed using the specific 
ionization energy loss and time-of-flight information, the ring-imaging Cherenkov technique and the 
kink-topology identification of weak decays of charged kaons. Transverse momentum spectra are 
measured from 0.1 up to 3 GeV/c for pions, from 0.2 up to 6 GeV/c for kaons and from 0.3 up 
to 6 GeV/c for protons. The measured spectra and particle ratios are compared with QCD-inspired 
models, tuned to reproduce also the earlier measurements performed at the LHC. Furthermore, the 
integrated particle yields and ratios as well as the average transverse momenta are compared with 
results at lower collision energies. 
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1 Introduction 

The majority of the particles produced at mid-rapidity in proton-proton collisions are low-momentum 
hadrons not originating from the fragmentation of partons produced in scattering processes with large 
momentum transfer. Their production, therefore, cannot be computed from first principles via pertur¬ 
bative Quantum Chromodynamics (pQCD). Currently available models describing hadron-hadron col¬ 
lisions at high energy, such as the event generators PYTHIA6 [1], PYTHIA8 [2, 3], EPOS [4, 5] and 
PHOJET [6], combine pQCD calculations for the description of hard processes with phenomenological 
models for the description of the soft component. The measurement of low-momentum particle produc¬ 
tion and species composition is therefore important as it provides crucial input for the modeling of the 
soft component and of the hadronization processes. Furthermore, it serves as a reference for the same 
measurement in Pb-Pb collisions to study the properties of the hot and dense strongly-interacting medium 
with partonic degrees of freedom, the Quark-Gluon Plasma (QGP), which is created in these collisions. 
In this paper, the measurement of primary n + ~, K + . p and p production at mid-rapidity in proton-proton 
collisions at ^/s = 7 TeV using the ALICE detector [7-10] is presented. Primary particles are defined 
as prompt particles produced in the collision including decay products, except those from weak decays 
of light flavour hadrons and muons. Pions, kaons and protons are identified over a wide momentum 
range by combining the information extracted from the specific ionization energy-loss (AE /Ax) mea¬ 
sured in the Inner Tracking System (ITS) [11] and in the Time Projection Chamber (TPC) [12], the time 
of flight measured in the Time-Of-Flight (TOF) detector [13], the Cherenkov radiation measured in the 
High-Momentum Particle Identification Detector (HMPID) [14] and the kink-topology identification of 
the weak decays of charged kaons. Similar measurements in proton-proton collisions at y/s = 900 GeV 
and 2.76 TeV are reported in [15-17] and are included, together with lower energy data [18-24], in the 
discussion of the evolution of particle production with collision energy. Similar measurements at the 
LHC have also been performed in the forward region [25]. 

The paper is organized as follows. In Section 2 the ALICE experimental setup is described, focusing 
on the detectors and the corresponding Particle IDentification (PID) techniques relevant for the present 
measurement. Details on the event and track selection criteria and the corrections applied to the measured 
raw yields are also presented. In Section 3 the results on the production of primary n ± , K + p and p are 
shown. These include the transverse momentum (pj) distributions and the pi-integrated production 
yields of each particle species and the K In and p In ratios. The evolution with collision energy of the p-\ - 
integrated particle yields, of their ratios and of their average transverse momenta (pj) is also presented. 
In Section 4 particle spectra and then - ratios (K In and pin) are compared with models, in particular with 
different PYTHIA tunes [1-3, 26, 27], EPOS [4, 5] and PHOJET [6]. Section 5 concludes the paper 
summarizing the results. 

2 Experimental setup and data analysis 

2.1 The ALICE detector 

The ALICE detector was specifically optimized to reconstruct and identify particles over a wide momen¬ 
tum range thanks to the low material budget, the moderate magnetic field and the presence of detectors 
exploiting ah the known PID techniques. A comprehensive description of the ALICE experimental setup 
and performance can be found in [7-10]. In the following, the PID detectors relevant for the analysis 
presented in this paper are briefly described, namely ITS, TPC, TOF and HMPID. They are located in the 
ALICE central barrel in a B = 0.5 T solenoidal magnetic field directed along the beam axis. The ITS, TPC 
and TOF detectors cover the full azimuth (<p) and have a pseudorapidity coverage of | t] < 0.9, while the 
HMPID covers the pseudorapidity interval |rj| < 0.55 and the azimuthal angle range 1.2° < <p < 58.5°. 
The ITS [11] is the innermost central barrel detector. It is composed of six cylindrical layers of silicon de¬ 
tectors, located at radial distances between 3.9 and 43 cm from the beam axis. The two innermost layers 
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are equipped with Silicon Pixel Detectors (SPD), the two intermediate ones are Silicon Drift Detectors 
(SDD), while the two outermost ones are Silicon Strip Detectors (SSD). The ITS provides high resolu¬ 
tion tracking points close to the beam line, which allows us to reconstruct primary and secondary vertices 
with high precision, to measure with excellent resolution the Distance of Closest Approach (DCA) of a 
track to the primary vertex, and to improve the track pj resolution. It is also used as a stand-alone tracker 
to reconstruct particles that do not reach the TPC or do not cross its sensitive areas. The SDD and SSD 
are equipped with analogue readout enabling PID via d E /dr measurements with a relative resolution of 
about 10 %. 

The TPC [12] is the main tracking detector of the ALICE central barrel. It is a large volume cylindri¬ 
cal chamber with high-granularity readout that surrounds the ITS covering the region 85 < r < 247 cm 
and -250 < z < +250 cm in the radial r and longitudinal z directions, respectively. It provides three- 
dimensional space points and specific ionization energy loss dE/dx with up to 159 samples per track. 
The relative dE/dx resolution is measured to be about 5.5% for tracks that cross from the centre of the 
outer edge of the detector. 

The TOF detector [13] is a large-area array of Multigap Resistive Plate Chambers (MRPC) with an in¬ 
trinsic time resolution of 50 ps, including the electronic readout contribution. It is a cylindrical detector 
located at a radial distance 370 < r < 399 cm from the beam axis. Particles are identified using simulta¬ 
neously the time-of-flight information with the momentum and track length measured with the ITS and 
the TPC. 

The HMPID [14] is a single-arm proximity-focusing Ring Imaging CHerenkov (RICH) detector located 
at 475 cm from the beam axis. The Cherenkov radiator is a 15-mm-thick layer of liquid C 6 F 14 (perfluo- 
rohexane) with a refractive index of n = 1.2989 at a photon wave length A = 175 nm, corresponding to 
a minimum particle velocity /3 m ; n = 0.77. 

In addition to the detectors described above that provide PID information, the VZERO system [28] is 
used for trigger and event selection. It is composed of two scintillator arrays, which cover the pseudora¬ 
pidity ranges 2.8 < 77 < 5.1 and —3.7 < 77 < —1.7. 

2.2 Data sample, event and track selection 

The results presented in this paper are obtained combining five independent analyses, namely ITS stand¬ 
alone, TPC-TOF, TOF, HMPID, kink, using different PID methods. The analyzed data are proton-proton 
collisions at ,/s = 7 TeV collected in 2010. During that period, the instantaneous luminosity at the 
ALICE interaction point was kept within the range 0.6-1.2x 10 29 cm 2 s~ 1 to limit the collision pile-up 
probability. Only runs with a collision pile-up probability smaller than 4% are used in this analysis, 
leading to an average pile-up rate of 2.5%. The number of events used in the five independent analyses 
is reported in Table 1 . The data were collected using a minimum-bias trigger, which required a hit in the 
SPD or in at least one of the VZERO scintillator arrays in coincidence with the arrival of proton bunches 
from both directions. This trigger selection essentially corresponds to the requirement of having at least 
one charged particle in 8 units of pseudorapidity. The contamination due to beam-induced background is 
removed off-line by using the timing information from the VZERO detector, which measures the event 
time with a resolution of about 1 ns, and the correlation between the number of clusters and track seg¬ 
ments (tracklets) in the SPD [15]. 

Selected events are further required to have a reconstructed primary vertex. For 87% of the triggered 
events, the interaction vertex position is determined from the tracks reconstructed in TPC and ITS. For 
events that do not have a vertex reconstructed from tracks, which are essentially collisions with low mul¬ 
tiplicity of charged particles, the primary vertex is reconstructed from the SPD tracklets, which are track 
segments built from pairs of hits in the two innermost layers of the ITS. Overall, the fraction of events 
with reconstructed primary vertex, either from tracks or from SPD tracklets, is of 91%. Accepted events 
are required to have the reconstructed vertex position along the beam direction, z, within ±10 cm from 
the centre of the ALICE central barrel. This ensures good rapidity coverage, uniformity of the particle 
reconstruction efficiency in ITS and TPC and reduction of the remaining beam-gas contamination. In the 
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Analysis 

# of events 

n 

K 

P 

ITS-sa 

5.4 10 7 

0.1-0.6 

0.2-0.5 

0.3-0.6 

TPC-TOF 

5.4-10 7 

0.25-1.2 

0.3-1.2 

0.45-2.0 

TOF 

5.4-10 7 

0.5-2.5 

0.5-2.4 

0.8-4.0 

HMPID 

OO 

© 

1.5-3.0 

1.5-3.0 

1.5-6.0 

Kink 

16.9 -10 7 

- 

0.2-6.0 

- 


Table 1: Number of analyzed events and px range (GeV/c) covered by each analysis. 

following analyses two different sets of tracks are used: the global tracks, reconstructed using informa¬ 
tion from both ITS and TPC, and the ITS-sa tracks, reconstructed by using only the hits in the ITS. To 
limit the contamination due to secondary particles and tracks with wrongly associated hits and to ensure 
high tracking efficiency, tracks are selected according to the following criteria. The global tracks are 
required to cross over at least 70 TPC readout rows with a value of X 2 /N clusters of the momentum fit in 
the TPC lower than 4, to have at least two clusters reconstructed in the ITS out of which at least one is 
in the SPD layers and to have a DCA to the interaction vertex in the longitudinal plane, DC A- < 2 cm. 
Furthermore, the daughter tracks of reconstructed kinks are rejected. This last cut is not applied in the 
kink analysis where a further px-dependent selection on the DCA of the selected tracks to the primary 
vertex in the transverse plane (DCA XV ) is requested. The global tracks that satisfy these selection criteria 
have a px resolution of 1% at p-\ = I GeV/c and 2% at p-\ = 10 GeV/c. The ITS-sa tracks are required to 
have at least four ITS clusters out of which at least one in the SPD layers and three in the SSD and SDD, 
X /^clusters < 2.5 and a DCA rv satisfying a /^-dependent upper cut corresponding to 7 times the DCA 
resolution. The selected ITS-sa tracks have a maximum px resolution of 6% for pions, 8% for kaons 
and 10% for protons in the px range used in the analysis. Both global and ITS-sa tracks have similar 
resolution in the DCA VV parameter that is 75 pm at px = 1 GeV/c and 20 pm at px = 15 GeV/c [29], 
which is well reproduced in the simulation of the detector performance. The final spectra are calculated 
for |y| < 0.5. 

2.3 Particle identification strategy 

To measure the production of 7i ± , K , p and p over a wide px range, results from five independent 
analyses, namely ITS-sa, TPC-TOF, TOF, HMPID and kink, are combined. Each analysis uses different 
PID signals in order to identify particles in the complementary px ranges reported in Table 1. In the 
following, the PID strategies used by ITS-sa, TPC-TOF and TOF analyses are briefly summarized since 
they are already discussed in detail in [15, 30], while the HMPID analysis, presented here for the first 
time, and the kink analysis, modified with respect to that described in [15], are presented in more detail. 

2.3.1 ITS stand-alone analysis 

In this analysis ITS-sa tracks are used and particles are identified by comparing the dE/dx measurement 
provided by the ITS detector with the expected values at a given momentum p under the corresponding 
mass hypotheses. In Fig. 1, the measured dE/dx values arc shown as a function of track momentum 
together with the curves of the energy loss for the different particle species, which are calculated using 
the PHOBOS parametrization [31] of the Bethe-Bloch curves at large /J y and with a polynomial to correct 
for instrumental effects. A single identity is assigned to each track according to the mass hypothesis for 
which the expected specific energy-loss value is the closest to the measured dE/dx for a track with 
momentum p. No explicit selection on the difference between the measured and expected values is 
applied except for a lower limit on pions set to two times the dE/dx resolution (a) and an upper limit on 
protons given by the mid-point between the proton and the deuteron expected dE/dx. The ITS dE/dx is 
calculated as a truncated mean of three or four dE/dx values provided by the SDD and SSD layers. The 
truncated mean is the average of the lowest two dE/dx values in case signals in all the four layers are 
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Fig. 1 : Distribution of AE/Ax as a function of momentum (p) measured in the ITS using ITS-sa tracks in |p| < 0.9. 
The continuous curves represent the parametrization of AE /dv for e, n, K and p while the dashed curves are the 
bands used in the PID procedure. 

available, or as a weighted average of the lowest (weight 1) and the second lowest (weight 1/2) values in 
the case where only three AE/Ax samples are measured. Even with this truncated mean approach, used to 
reduce the effect of the tail of the Landau distribution at large AE / Ax, the small number of samples results 
in residual non-Gaussian tails in the AE/ dr distribution, which are partially reproduced in simulation. 
These non-Gaussian tails increase the misidentification rate, e.g. pions falling in the kaon identification 
bands. The misidentification probability is estimated using a Monte-Carlo simulations where the particle 
abundances were adjusted to those observed in data. This correction is at most 10% in the pj range 
of this analysis. In order to check possible systematic effects due to these non-Gaussian tails and their 
imperfect description in Monte-Carlo simulations, the analysis was repeated with different strategies for 
the particle identification, namely using a 3a compatibility band around the expected AE/Ax curves and 
extracting the yields of pions, kaons and protons using the unfolding method described in [15], which is 
based on fits to the AE/Ax distributions in each pj interval. The difference among the results from these 
different analysis strategies is assigned as a systematic uncertainty due to the PID. 

2.3.2 TPC-TOF analysis 

In this analysis global tracks are used and particle identification is performed by comparing the measured 
PID signals in the TPC and TOF detectors (AE / Ax, time of flight) with the expected values for different 
mass hypotheses. An identity is assigned to a track if the measured signal differs from the expected 
value by less than three times its resolution a. For pions and protons with pj < 0.6 GeV/c and kaons 
with pj < 0.5 GeV/c, a compatibility within 3a is required on the AE/Ax measurement provided by 
the TPC computed as a truncated mean of the lowest 60% of the available AE/Ax samples. The AE/Ax 
resulting from this truncated mean approach is Gaussian and it is shown in Fig. 2 as a function of the 
track momentum together with the expected energy-loss curves (see [32] for a discussion about the 
AE/Ax parametrization). Above these pj thresholds, i.e. pj > 0.6 GeV/c for pions and protons and 
pj > 0.5 GeV/c for kaons, a three a requirement is applied to both the AE/Ax measurement provided 
by the TPC and the time of flight f to f provided by the TOF detector. The time of flight t to f> as will be 
described in more detail in the next section, is the difference between the arrival time Ttof measured with 
the TOF detector and the event start time to, namely t to f = Ttof — to- The additional condition on the TOF 
signal helps in extending the particle identification on a track-by-track basis to higher pj where the TPC 
separation power decreases. The particles for which the time-of-flight signal is available are a sub-sample 
of the global tracks reconstructed using ITS and TPC information. The time-of-flight information is not 
available for tracks that cross inactive regions of the TOF detector, for particles that decay or interact 
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Fig. 2: Distribution of dE/dv as a function of momentum (p) measured in the TPC using global tracks for \rj\ <0.9. 
The continuous curves represent the Bethe-Bloch parametrization. 




Fig. 3: Particle velocity /3 measured by the TOF detector as a function of the rigidity p/z, where z is the particle 
charge, for 1 77 1 < 0.9. 

with the material before the TOF and for tracks whose trajectory, after prolongation from the TPC outer 
radius, is not matched with a hit in the TOF detector. The fraction of global tracks with associated 
time-of-flight information (TOF matching efficiency) depends on the particle species and pj as well as 
on the fraction of the TOF active readout channels. For the data analysis presented in this paper the 
matching efficiency increases with increasing pj until it saturates, e.g. at about 65% for pions with pj 
> 1 GeV/c. In Fig. 3 the velocity /J of the tracks, computed from the trajectory length measured with 
the ITS and TPC and the time of flight measured with the TOF, is reported as a function of the rigidity 
p/z, where z is the charge assigned based on the measured direction of the track curvature. More than 
one identity can be assigned to a track if it fulfils PID and rapidity selection criteria for different particle 
species. The frequency of such cases is at most 0.5% in the momentum range used in this analysis. The 
misidentification of primary particles is computed and corrected for using Monte-Carlo simulations. It 
is at most 2% for pions and protons and 8 % for kaons in the considered pj ranges. The correction of the 
raw spectra for the misidentified particles provides also a way to remove the overestimation of the total 
number of particles introduced by the possibility, described above, to assign more than one identity to a 
track. 
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Fig. 4: Distribution of At,- assuming the pion mass hypothesis in the transverse momentum interval 
1.9 < p'Y < 2.0 GeV/c. The data (black points) are fitted with a function (light blue line) that is the sum of 
the signal due to pions (green dotted line) and the two populations corresponding to kaons (red dotted line) and 
protons (purple dashed line). 

2.3.3 TOF analysis 

This analysis uses the sub-sample of global tracks for which a time-of-flight measurement is available. 
The PID procedure utilizes a statistical unfolding approach that provides a pj reach higher than the three 
a approach described in the previous section. The procedure is based on the comparison between the 
measured time of flight from the primary vertex to the TOF detector, t to f, and the time expected under 
a given mass hypothesis, ? ; exp (i = n, K, p), namely on the variable A?,• = t to f— ff xp . As mentioned in 
the previous section, the time of flight f to f is defined as the difference between the time measured with 
the TOF detector Ttof and the event start time to- The to value is computed from the analyzed tracks 
themselves on an event-by-event basis, using a combinatorial algorithm which compares the measured 
Ttof with the expected ones for different mass hypotheses. The track under study is excluded to avoid 
any bias in the PID procedure [13, 15]. In case the TOF to algorithm fails, the average beam-beam 
interaction time is used. The former approach provides a better to resolution, but requires at least three 
reconstructed tracks with an associated TOF timing measurement. The yield of particles of species i in 
a given pj interval is obtained by fitting the distribution of the variable At, obtained from all the tracks 
regardless of the method used to compute the to- This distribution is composed of the signal from particles 
of species i, which is centered at At,- = 0, and two distinct populations corresponding to the other two 
hadron species, j,k / i. The At,- distribution is therefore fitted with the sum of three functions /(At,), one 
for the signal and two for the other hadron species, as shown in Fig. 4. The /(At,) functional forms are 
defined using the data in the region of clear species separation. The TOF signal is not purely Gaussian 
and it is described by a function /(At,) that is composed of a Gaussian term and an exponential tail at 
high At, mainly due to tracks inducing signals in more than one elementary detector readout element [13]. 
The raw yield of the species i is given by the integral of the signal fit function. The reach in p-\ of this PID 
method depends on the resolution of At, that is the combination of the TOF detector intrinsic resolution, 
the uncertainty on the start time and the tracking and momentum resolution. Its value, for the data used 
in this analysis, is about 120 ps leading to 2a pion-kaon and kaon-proton separation at p-\ = 2.5 GeV/c 
and pj = 4.0 GeV/c, respectively. This PID procedure has the advantage of not requiring a Monte-Carlo 
based correction for misidentification because the contamination under the signal of particles of species 
i due to other particle species is accounted for by the background fit functions. 
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Fig. 5: Display of a Cherenkov ring detected in a module of HMPID for an inclined track crossing the detector. 
The colors are proportional to the pad charge signal. 
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2.3.4 HMPID analysis 

The HMPID is a Ring Imaging Cherenkov (RICH) detector in a proximity focusing layout in which 
the primary ionizing charged particle generates Cherenkov light inside a liquid CgFi 4 radiator [14]. 
The UV photons are converted into photoelectrons in a thin Csl him of the PhotoCathodes (PCs) and 
the photoelectrons are amplified in an avalanche process inside a Multi-Wire Proportional Chamber 
(MWPC) operated with CH 4 . To obtain position sensitivity for the reconstruction of the Cherenkov 
rings, the PCs are segmented into pads. The final image of a Cherenkov ring is then formed by a cluster 
of pads (called “MIP” cluster) associated to the primary ionization of the particle and the photoelectron 
clusters associated to Cherenkov photons. In Fig. 5 a typical Cherenkov ring is shown. 

In this analysis, the sub-sample of global tracks that reach the HMPID detector and produce the Cherenkov 
rings is used. Starting from the photoelectron cluster coordinates on the photocathode, a back-tracking 
algorithm calculates the corresponding single photon Cherenkov angle by using the impact angle of a 
track extrapolated from the central tracking detectors up to the radiator volume. A selection on the dis¬ 
tance (r/MiP-trk) computed on the cathode plane between the centroid of the MIP cluster and the track 
extrapolation, set to r/MiP-trk < 5 cm, rejects fake associations in the detector. Background discrimina¬ 
tion is performed using the Hough Transform Method (HTM) [33]. 

The mean Cherenkov angle (P c kov) is obtained if at least three photoelectron clusters are detected. For 
a given track, (P c kov) is computed as the weighted average of the single photon angles (if any) selected 
by HTM. Pions, kaons and protons become indistinguishable at high momentum when the resolution on 
(0 c kov) reaches 3.5 mead. The angle (0 c kov) as a function of the track momentum is shown in Fig. 6 , 
where the solid lines represent the 0 c k O v dependence on the particle momentum 


dckm — COS 


s/p 2 + m 2 
np 


( 1 ) 


where n is the refractive index of the liquid radiator, m the mass of the particle and p its momentum. 
This analysis is performed for p >1.5 GeV/c, where pions, kaons and protons produce a ring with enough 
photoelectron clusters to be reconstructed. If the track momentum is below the threshold to produce 
Cherenkov photons, background clusters could be wrongly associated to the track. As an example the few 
entries visible in Fig. 6 between the pion and kaon bands at low (P c kov) correspond to wrong associations 
of clusters with a kaon or a proton below the threshold to produce Cherenkov photons. 

The particle yields are extracted from a ht to the Cherenkov angle distribution in narrow transverse 
momentum intervals. In Fig. 7, examples of the reconstructed Cherenkov angle distributions in two 
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Fig. 6: Mean Cherenkov angle (0 c kov) measured with HMPID in its full geometrical acceptance as a function of 
the particle momentum p for positively and negatively charged tracks. The solid lines represent the theoretical 
curves for each particle species. 

narrow p \ intervals (3.4 < pj < 3.6 GeV/c and 5 < p-\ < 5.5 GeV/c) for negatively-charged tracks are 
shown. 

The background, mainly due to noisy pads and photoelectron clusters from other rings overlapping to 
the reconstructed one, is negligible in the momentum range considered in this analysis. The fit function 
(shown as a solid line in Fig. 7) is a sum of three Gaussian functions, one for each particle species (dashed 
lines), whose mean and sigma are fixed to the Monte-Carlo values. The extracted separation power of 
hadron identification in the HMPID as a function of pj is shown in Fig. 8. The separation between 
pions and kaons (kaons and protons) is expressed as the difference between the means of the (0 c kov) 
angle Gaussian distributions for the two given particle species (A^k or Ak, p ) divided by the average of 
the Gaussian widths of the two distributions, i.e. (a^+a k)/2 or (aK+a p )/2. A separation at 3a level in 
(0 c kov) is achieved up to p-\ = 3 GeV/c for K -n and up to p T = 5 GeV/c for K-p. The separation at 
6 GeV/c for K-p can be extrapolated from the curve and it is about 2.5a. The HMPID geometrical 
acceptance is about 5% for tracks with high momentum. Therefore the analysis of HMPID required to 
analyze a larger data sample with respect to the other PID methods, as reported in Table 1. The total 
efficiency is the convolution of the tracking, matching and PID efficiencies. The PID efficiency of this 
method is determined by the Cherenkov angle reconstruction efficiency. It has been computed by means 
of Monte-Carlo simulations and it reaches 90% for particles with velocity /3 ~ 1. As a cross check, the 
PID efficiency has been determined using clean samples of protons and pions from A and K° decays. 
The measured efficiency agrees within the statistical uncertainties with the Monte-Carlo estimates, in the 
momentum range 1.5 < p-\ < 6 GeV/c. Moreover, the correction due to the rfMiP-trk cut is computed 
from the same sample of identified protons and pions from A and K ( s ' decays. 

2.3.5 Kink analysis 

Charged kaons can also be identified in the TPC by reconstructing their - weak-decay vertices, which 
exhibit a characteristic kink topology defined by a decay vertex with two tracks (mother and daughter) 
having the same charge. This procedure extends the measurement of charged kaons on a track-by-track 
basis to pj = 6 GeV/c. The algorithm for the kink reconstruction is applied inside a fiducial volume of 
the TPC, namely 130 < R < 200 cm, needed to reconstruct both the mother and the daughter tracks. The 
mother track is selected with similar criteria as the global tracks (Section 2.2), but with a looser selection 
on the minimum number of TPC clusters, which is set to 20, and a wider rapidity range set to |y| < 0.7 to 
increase the statistics of kink candidates. No selections are applied on the charged daughter track. The 
reconstructed invariant mass M u v is calculated assuming the charged daughter track to be a muon and the 
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Cherenkov angle <# ckov ) (rad) 

Fig. 7: Distributions of (0 c kov) measured with the HMPID in the two narrow p-\- intervals 3.4 < p-\ < 3.6 GeV/c 
(top) and 5 < pt < 5.5 GeV/c (bottom) for tracks from negatively-charged particles. Solid lines represent the total 
fit (sum of three Gaussian functions). Dotted lines correspond to pion, kaon and proton signals. The background 
is negligible. 



Fig. 8: Separation power (n a ) of hadron identification in the HMPID as a function of p-\ . The separation of 
pions and kaons (kaons and protons) is defined as the difference between the average of the Gaussian distributions 
of (0 c kov) f°r the two hadron species divided by the average of the Gaussian widths of the two distributions. 
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Fig. 9: Kink invariant mass M^ v in data (red circles) and Monte-Carlo (black line) for summed particles and 
antiparticles, integrated over the mother transverse momentum range 0.2 < pj < 6.0 GeV/c and |v < 0.7 before 
(top panel) and after (bottom panel) the topological selections, based mainly on the qj and the maximum decay 
opening angle. 

undetected neutral daughter track to be a neutrino. The neutrino momentum is the difference between the 
measured momenta of the mother particle and of the charged daughter. The M^ v distribution, for summed 
positive and negative charges, integrated over the mother transverse momentum range 0.2 < pj < 6.0 
GeV/c is reported in the top panel of Fig. 9 for both data and PYTHIA simulations normalized to the 
same number of entries. Three peaks are present: one centered on the kaon mass due to the kaon decays 
K—» p + Vu (branching ratio BR = 63.55%), one centered at M^ v = 0.43 GeV/c 2 due to the K —> n + n° 
decay (BR = 20.66%), whose kinematics is calculated with wrong mass assumptions, and the peak due 
to pion decays n —> p -f (BR = 99.99%). The width of the peaks reflects the momentum resolution 
of the detector, which is well reproduced in Monte-Carlo simulations. The two-body kinematics of the 
kink topology allows one to separate kaon decays from the main source of background due to charged 
pion decays [15]. In the p + v^ channel, the upper limit of the q-\- variable, where q-\- is defined as 
the transverse momentum of the daughter track with respect to the mother’s direction, is 236 MeV/c for 
muons from kaon decays and 30 MeV/c for muons from pion decays. To remove most of the pion decays, 
a qj > 120 MeV/c selection is applied. The background is further reduced by rejecting kink decays for 
which the decay angle, namely the angle between the momenta of the mother and the charged daughter 
tracks is larger than the maximum angle allowed in the hypothesis K— > p + v f ,. The bottom panel of 
Fig. 9 shows the invariant mass distribution of the kaon candidates with mother transverse momentum 
0.2 < pt < 6.0 GeV/c after the topological selection criteria for kaon identification (mainly the q\ and 
decay angle cuts) are applied. It is evident that only the two peaks coming from kaon decays are present 
while the pion background peak is removed. The broad structure on the left originates from the three- 
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body decays of kaons. The agreement between data and simulations in this figure (Fig. 9) is better than 
8 %. Most of the selected mother tracks have ad E/dx in the TPC which is compatible with the values 
expected for kaons. Tracks outside 3.5a from the expected kaon d E/dx have been removed to attain a 
purity > 97 % in the pj range studied in this analysis. These rejected tracks are less than 4%, have pj < 
0.8 GeV/c and are, according to Monte-Carlo studies, pions. The raw kaon spectra are obtained from the 
integral of the invariant mass distribution computed in narrow p-\ intervals after the topological selection 
criteria on the qj , the decay opening angle, and the compatibility with the expected dE/dx for kaons are 
applied. The kaon misidentification is computed and corrected for by using Monte-Carlo simulations. 
It depends on the mother’s transverse momentum with a maximum value of 3.6% at 0.8 GeV/c and a 
minimum of 2% at 1 GeV/c, remaining almost flat up to pj = 6 GeV/c. Its average value in the pj range 
considered in this analysis is 2.1%. 

2.4 Correction of raw spectra 

To obtain the pj distributions of primary n, K and p, the contribution of secondaries is subtracted from 
the raw spectra. Then, the spectra are corrected for the PID efficiency, the misidentification probability, 
the acceptance, the reconstruction and the selection efficiencies according to 


d 2 N 
dpjdy 


— N ]raw (PT ) 


1 

ApjAy 


1 — s { pt ) 
£(pt) 


■Apt), 


( 2 ) 


where N Taw (pj) Aj> |. Av is the raw yield in a given pj interval, s(pj) is the total contamination including 
effects of secondary and misidentified particles, £(pj) is the acceptancexefficiency including PID effi¬ 
ciency, detector acceptance, reconstruction and selection efficiencies and f(pj) is an additional factor to 
correct for imperfections of the cross sections for antiparticle interactions with the material used in the 
GEANT3 code. 


The contamination due to weak decays of light flavour hadrons (mainly K* 1 affecting n spectra and A and 
affecting p spectra) and interactions with the material has to be computed and subtracted from the 
raw spectra. Since strangeness production is underestimated in the event generators and the interactions 
of low pj particles with the material are not properly modelled in the transport codes, the secondary- 
particle contribution is evaluated with a data-driven approach. This approach exploits the high resolution 
determination of the track impact parameter in the transverse plane, DCA XV , and the fact that secondary 
particles from strange hadron decays and interactions with the detector material, originate from sec¬ 
ondary vertices significantly displaced from the interaction point and, therefore, their tracks have, on 
average, larger absolute values of DCA X> , with respect to primary particles. Hence, for each of the PID 
techniques described in the previous sections, the contribution of secondary particles to the measured raw 
yield of a given hadron species in a given pj interval is extracted by fitting the measured distributions of 
DCAjrv of the tracks identified as particles of the considered hadron species. The DCA xy distributions are 
modelled with three contributions, called templates. Their shapes are extracted for each pj interval and 
particle species from simulations, as described in [30], and represent the DCA a v distributions of primary 
particles, secondary particles from weak decays of strange hadrons and secondary particles produced in 
interactions with the detector material, respectively. An example for protons in the interval 0.55 < pj < 
0.60 GeV/c is shown in Fig. 10. The correction for secondary particle contamination is relevant for n + 
(from 10% at low pj to less than 2% at high pj), p and p (from 35% at low pj to 2% at high pj). Due 
to the different track and PID selections the contribution of secondaries is different for each analysis. 

In the case of kaons, the contamination from secondary particles is negligible, except for the TPC- 
TOF analysis where a contamination originating from secondary e* produced by photon conversions 
in the detector material is present. This contamination is significant only in the momentum range 
0.4 < p < 0.6 GeV/c, where the dE/dx of kaons and electrons in the TPC gas are similar, not al¬ 
lowing for their separation, as shown in Fig. 2. Therefore, in the case of kaons, the fit to the DCA rv 
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Fig. 10: Proton DCA rv distribution in the range 0.55 < p \ < 0.60 GeV/c together with the Monte-Carlo templates 
for primary protons (green dotted line), secondary protons from weak decays (red dotted line) and secondary 
protons produced in interactions with the detector material (blue dashed line) which are fitted to the data. The light 
blue line represents the combined fit, while the black dots are the data. 
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Fig. 11: Correction factors (e(pj) in Eq. 2) for K ', K + and p (left panel) and their antiparticles (right panel) 
accounting for PID efficiency, detector acceptance, reconstruction and selection efficiencies for ITS-sa (red circles), 
TPC-TOF (light blue squares), TOF (green diamonds), HMPID (black stars) and kink (purple crosses) analyses. 
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distributions is used only in the TPC-TOF analysis for pt <0.5 GeV/c to subtract the contamination due 
to secondary e + . This contamination is about 16% for p-\ = 0.5 GeV/c. 

The resulting spectra are corrected for the detector acceptance and for the reconstruction and selection 
efficiencies. This correction is specific to each analysis and accounts for the acceptance of the detector 
used in the PID procedure, the trigger selection and the vertex and track reconstruction efficiencies. 
They are evaluated by performing the same analyses on simulated events generated with PYTHIA 6.4 
(PerugiaO tune) [26]. The particles are propagated through the detector using the GEANT3 transport 
code [34], where the detector geometry and response as well as the data taking conditions are reproduced 
in detail. 

In Fig. 11 the efficiency e(pj), specific to each analysis, accounting for PID efficiency, acceptance, 
reconstruction and selection efficiencies are shown. The lower value of e for HMPID and kink analyses 
is due to the limited geometrical acceptance of the HMPID detector and to the limited TPC fiducial 
volume used for the kink vertex reconstruction. The drop in the correction for the TPC-TOF analysis at 
Pt = 0.6 GeV/c for pions and protons and p-\ =0.5 GeV/c for kaons is due to the efficiency of track 
propagation to the TOF. The ITS-sa analysis has a larger kaon efficiency than TPC-TOF analysis at low 
Pt because the ITS-sa tracking allows the reconstruction of kaons that decay before reaching the TPC. 
The corrections for particles (left panel of Fig. 11) and antiparticles (right panel) are compatible within 
the uncertanties. 

Since GEANT3 does not describe well the interaction of low-momentum p and K with the material, 
corrections to the efficiencies, estimated with a dedicated FLUKA simulation [30, 35], are applied. The 
correction factor f(pj) is 0.71 < f(pj) < 1 for p and 0.95 < f(pj) < 1 for K . 

The corrected spectra are, finally, normalized to the number of inelastic proton-proton collisions that is 
obtained from the number of analyzed minimum-bias events via the scaling factor 0.852 as described 
in [36]. 

2.5 Systematic uncertainties 

The main sources of systematic uncertainties, for each analysis, are summarized in Table 2. They are 
related to the PID procedure, the subtraction of the contribution from secondary particles, imperfect 
description of the material budget in the Monte-Carlo simulation, particle interactions in the detector 
material, tracking efficiency and variables used for the track selection. 

The systematic uncertainties introduced by the PID procedure are estimated differently depending on 
the specific analysis. In the ITS-sa analysis different techniques are used for the identification: a 3a 
compatibility cut and an unfolding method as described in section 2.3.1. In the TPC-TOF analysis the 3a 
selection is varied to 2a and 4a. Furthermore, the systematic uncertainty on the estimated contamination 
from misidentified hadrons, which is due to the different relative abundances of pions, kaons and protons 
in data and simulation, has been estimated to be below 1% for pions and protons and below 4% for kaons. 
In TOF and HMPID analyses the parameters of the fit function used to extract the raw yields are varied 
(one at a time) by ± 10%. 

The systematic uncertainty due to the subtraction of secondary particles is estimated by changing the fit 
range of the DCA av distribution. The shape of the DCA X) , template for p and p from weak decays is also 
varied by modifying the relative contribution of the different mother particles. The main sources of p 
and p from weak decays are A and £ + (and their antiparticles), which have significantly different mean 
proper decay lengths (ct = 7.89 cm and 2.404 cm respectively [37]). Therefore, the DCA template of 
protons from weak decays depends on the A to £ + ratio in the event generator used in the simulation. 

To evaluate the systematic effect due to the uncertainty in the material budget (about ±7% [38]), the 
efficiency corrections are computed by using Monte-Carlo simulations with the material budget modified 
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Source of uncertainty common to all the analyses 

7l ± 

K± 

p and p 


Correction for secondaries 

Material budget 

Cross sections for interactions in the material 

ITS-TPC matching (excluded in ITS-sa analysis) 

< 1% 

5%-Negl. 

2% — 1% 

3%- 

3%-Negl. 
4% — 1% 

3% 

5% — 1.5% (p) 
1.5% (p) 

3% - Negl. 

4% - Negl. (p) 

6% — 1% (p) 

3% 

Source of uncertainty specific to an analysis 

7l ± 

K ± 

p and p 

ITS-sa PID 

2% 

4% 

4.5% 

Tracking efficiency (ITS-sa tracks) 

3% 

3% 

3% 

ExB effect 

3% 

3% 

3% 

TPC-TOF PID 

< 1% 

1% — 5% 

< 1% 

Tracking efficiency (global tracks) 

2% 

2% 

2% 

Matching efficiency 

3% 

6% 

4% 

(px > 0.5 GeV/c for K and 0.6 GeV/c for n, p) 




TOF PID 

0.5% - 3% 

1% —11% 

1% —14% 

Tracking efficiency (global tracks) 

2% 

2% 

2% 

Matching efficiency 

3% 

6% 

4% 

HMPID PID 

4% 

5% 

5%-9% 

Tracking efficiency (global tracks) 

5% 

5% 

7% 

dMIP-trk cut 

2% - 6% 

2% - 6% 

2% - 6% 

Kink PID 


3% 


Tracking efficiency (global tracks) 


2% 


Kink reconstruction efficiency 


3% 


Kink contamination 


3.6%-2% 


Table 2: Sources of systematic uncertainties on the corrected spectra . In case of pj-dependent systematic 

uncertainty, the values in the lowest and highest p-\ intervals are reported. 

by this percentage. The systematic uncertainties in modelling the particle interactions with the detector 
material are evaluated using different transport codes, as described in [30]. 

For all the analyses, the systematic uncertainties related to tracking procedure are estimated by varying 
the track selection criteria (e.g. number of crossed readout rows in TPC, number of clusters in ITS, 


DCAj, DCA ty ) reported in Section 2.2. For global tracks an additional uncertainty, related to the ITS- 
TPC matching, is also included. It is estimated by comparing the matching efficiency in data and Monte- 
Carlo simulations. 

Further systematic uncertainty sources, specific to each analysis, are also evaluated. In the case of the 
ITS-sa analysis, the Lorentz force causes shifts of the cluster position in the ITS, pushing the charge in 
opposite directions depending on the polarity of the magnetic field of the experiment (E xB effect). This 
effect is not fully reproduced in the simulation. It is estimated by analysing data samples collected with 
different magnetic field polarities, which resulted in an uncertainty of 3%. In the case of TPC-TOF and 
TOF analyses, the influence of the material budget on the matching of global tracks with hits in the TOF 
detector is computed by comparing the matching efficiency for tracks traversing a different amount of 
material, in particular sectors with and without Transition Radiation Detector (TRD) modules installed. 
In the HMPID analysis, the ^MiP-trk cut selection is varied to check its systematic effect on the matching 
of global tracks with HMPID signals. 

In the kink analysis, the total systematic uncertainty is calculated as the quadratic sum of the contribu¬ 
tions listed in Table 2. The kaon misidentification correction (I - purity) described in Section 2.3.5, which 
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is on average 2.1%, depends on the relative particle abundances in the Monte-Carlo and a pr-dcpcndent 
uncertainty of about 2% on the purity is estimated. The kink identification uncertainty (3%, almost flat 
in the considered p-\ region) is also estimated with Monte-Carlo simulations by comparing the results by 
varying slightly some parameters of the analysis: the fiducial volume of the TPC is increased from the 
nominal 130 < R < 200 cm tol20</?<210 cm, the qj threshold is reduced from the nominal 120 
MeV/c to 40 MeV/c, and the requirement on the number of TPC clusters of the mother track is increased 
from the nominal 20 to 50 clusters. 

The systematic uncertainty on the efficiency for findable kink vertices was estimated to be 3% indepen¬ 
dently of pj by comparing, in real data and Monte Carlo simulations, the number of raw reconstructed 
kinks per kink radius unit in two different fiducial volumes inside the TPC, namely 130-200 cm and 
140-190 cm. 

Finally, a systematic uncertainty common to each analysis is related to the normalization to inelastic 
collisions. The normalization factor was evaluated in [36] and it is 0.852 ^qq|q. 

Ah described uncertainties are strongly correlated among the pj bins. Most of the uncertainties (e.g. 
tracking efficiency, ITS-TPC matching, TOF matching, material budget or PID) are also correlated be¬ 
tween the different particle species. 

3 Results 

The mid-rapidity (|y| < 0.5) transverse momentum spectra of n ] +n 4 K 1 +K and p+p obtained with 
the five analysis techniques discussed in Section 2, normalized to the number of inelastic collisions 
/V[\e[ , are reported in the top panel of Fig. 12. For a given hadron species, the spectra of particles and 
antiparticles are found to be compatible within uncertainties. Therefore, all the spectra shown in this 
section arc reported for summed charges. Since in their overlap p \ regions the spectra from the different 
PID techniques are consistent within uncertainties, they are averaged in a sequential procedure. The first 
step consists in averaging the two analyses whose results are the most closely correlated (namely TPC- 
TOF and TOF). Successively, the other analyses are added one-by-one to the running average according 
to then - degree of correlation with the previous ones. At each step of this sequential procedure, a weighted 
average of two spectra is computed by using as weights the inverse of the squares of the uncorrelated 
systematic uncertainties. The uncorrelated and correlated uncertainties are propagated separately through 
the weighted average formula. In Fig. 13 the n, K and p spectra, resulting from the combination of the 
five analyses, are reported. The bottom panels of Fig. 12 show the ratios between the spectra from 
each analysis and the combined one: the former are considered with their total systematic uncertainties 
while the latter without uncertainty. The uncertainty due to the normalization to inelastic collisions 
(^ 4 %), common to the five PID analyses, is not included. The agreement between each analysis and the 
combined one is satisfactory, being within the total systematic uncertainties. 

To extrapolate to zero and infinite momentum, the combined spectra reported in Fig. 13 are fitted with 
the Levy-Tsallis function [39, 40] 


d 2 N 
dpxdy 


d N , 

PT -* (1 + 


mj — mo 
nC 


(3) 


where 

K = (»-!)(»- 2) 
nC (nC + mo (n — 2)) 

mj = \JPj + wJq, //jo is the particle rest mass and C, n and the yield dA/dy are the free parameters. The 

Levy-Tsallis function describes rather well the spectra. The x 2 per number of degrees of freedom (ndf) 
of the fit are lower than unity (see Table 3) due to residual correlations in the point-to-point systematic 
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Fig. 12: Top panel: p-\ spectra of n, K and p, sum of particles and antiparticles, measured with ALICE at mid¬ 
rapidity (|y| < 0.5) in pp collisions at sfs = 7 TeV by using different PID techniques. The spectra are normalized 
to the number of inelastic collisions. Statistical (vertical error bars) and systematic (open boxes) uncertainties are 
reported. The horizontal width of the boxes represents the px-bin width. The markers are placed at the bin centre. 
Bottom panels: ratio between the spectra obtained from each analysis and the combined one. The error bands 
represent the total systematic uncertainties for each analysis. The uncertainty due to the normalization to inelastic 
collisions (+ 4 %), common to the five PID analyses, is not included. 

uncertainties. In Table 3 the values of the px-integrated yield d/V/dy and of the mean transverse mo¬ 
mentum (px) are reported for each particle species. They are obtained using the measured data in the pj 
range where they are available and the Levy-Tsallis function fitted to the data elsewhere, to extrapolate 
to zero and infinite momentum. The lowest px experimentally accessible and the fraction of yield con¬ 
tained in the extrapolated region are also reported in the table. The extrapolation to infinite momentum 
gives a negligible contribution to the values of both dN/dy and (pt). The d/V /dy and (px) uncertainties 
reported in Table 3 are the combination of the statistical and the systematic ones. The statistical uncer¬ 
tainties are negligible, while the systematic uncertainties are the sum of two independent contributions. 
The first contribution is due to the systematic uncertainties on the measured px-differential yields and 
it was estimated by repeating the Levy-Tsallis fits moving the measured points within their systematic 
uncertainties. The second contribution is due to the extrapolation to zero momentum and it is estimated 
using different fitting functions (namely modified Hagedorn [41] and UA1 parametrization [42]). Results 
for positively and negatively charged particles, separately, are also reported. It should be noticed that the 
yields of particles and antiparticles are compatible within uncertainties. 

In Fig. 14 the px spectra of identified charged hadrons, sum of particles and antiparticles, measured 
with ALICE at y/s = 7 TeV are compared to the results obtained by the CMS Collaboration at the same 
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Fig. 13: Combined p-\ spectra of n, K and p, sum of particles and antiparticles, measured with ALICE at mid¬ 
rapidity (|y| < 0.5) in pp collisions at y/s = 7 TeV normalized to the number of inelastic collisions. Statistical 
(vertical error bars) and systematic (open boxes) uncertainties are reported. The uncertainty due to the normaliza¬ 
tion to inelastic collisions (+ 4 %) is not shown. The spectra are fitted with Levy-Tsallis functions. 


Particle 

cLV//dy 

<Pt> 

(GeV/c) 

*7«df 

L .p T 
(GeV/c) 

Extr. 

% 

K + + 

4.49 ±0.20 

0.466 ±0.010 

19.1/38 

0.10 

9 

K+ + KL 

0.572 ±0.032 

0.773 ±0.016 

5.0/45 

0.20 

10 

p + p 

0.247 ±0.018 

0.900 ±0.029 

10.8/43 

0.30 

12 

n + 

2.26±0.10 

0.464 ±0.010 

24.0/38 

0.10 

9 

n~ 

2.23 ±0.10 

0.469 ±0.010 

15.0/38 

0.10 

9 

K+ 

0.286±0.016 

0.777 ±0.016 

7.4/45 

0.20 

9 

K± 

0.286±0.016 

0.770±0.016 

10.0/45 

0.20 

10 

P 

0.124 ±0.009 

0.900 ±0.027 

9.5/43 

0.30 

12 

P 

0.123 ±0.010 

0.900 ±0.032 

12.3/43 

0.30 

12 


Table 3: d/V/dy and (p-y) extracted from Levy-Tsallis fits to the measured n, K, p spectra in inelastic pp collisions 
at yj~s = 7 TeV for |y| < 0.5 with combined statistical and systematic uncertainties (statistical uncertainties are 
negligible) together with the pj of the lowest experimentally accessible point (L. p T ) and the extrapolated fraction. 
The systematic uncertainty on d/V/dv due to normalization to inelastic collisions (_ 4 %) is not included. 


centre-of-mass energy [17]. Even though the measurements are performed in different rapidity intervals 
(|y| < 0.5 for ALICE, jy < 1 for CMS), they can be compared since the pj spectra are essentially 
independent of rapidity for |y| < 1. A similar comparison at y/s = 0.9 TeV is reported in [17]. At both 
energies, the ALICE spectra are normalized to the number of inelastic collisions, while the CMS results 
are normalized to the double-sided selection (at least one particle with E > 3 GeV in both —5 < 77 < —3 
and 3 < rj < 5). An empirical scaling factor of 0.78, computed by the CMS Collaboration in [17] for 
the spectra measured in pp collisions at y/s = 0.9 TeV, is therefore applied to the CMS data points at y/s 
= 7 TeV, to take into account the different event selections (details are given in [17]). With this scaling, 
the pion and kaon spectra measured with ALICE and CMS are found to agree within uncertainties. 
The proton spectra have different slopes: for pj < 1 GeV/c the ALICE and CMS results agree within 
uncertainties, while at higher pj a discrepancy of up to 20% is observed. 

In Fig. 15 the 7r, K and p integrated yields, d N /d y, are compared with similar measurements in the central 
rapidity region at various collision energies. In particular, results from ALICE at y/s = 900 GeV [15] 
and = 2.76 TeV [16], PHENIX at y/s = 62.4 GeV and = 200 GeV [18] and CMS, scaled by the 


18 












Pion, kaon and proton productionin pp collisions at \fs -1 TeV 


ALICE Collaboration 



Fig. 14: Comparison of pj spectra of n, K and p (sum of particles and antiparticles) measured by the ALICE 
(|y| < 0.5) and CMS Collaborations (|v| < 1) in pp collisions at v (v = 7 TeV. The CMS data points are scaled by 
the empirical factor 0.78, as described in [17]. Inset plot: Ratios between ALICE and CMS data in the common pj 
range. The combined ALICE and CMS statistical (vertical error bars) and systematic (open boxes) uncertainties 
are reported. The combined ALICE (+ 4 %) and CMS (±3%) normalization uncertainty is shown as a grey box 
around 1 and not included in the point-to-point uncertainties. 
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Fig. 15: px-integrated yields d/V/dy of n, K and p as a function of the centre-of-mass energy in pp collisions. 
PHENIX results are for separate charges while CMS and ALICE results are the average of the d/V/d v of particles 
and antiparticles. ALICE and CMS points are slightly shifted along the x-axis for a better visualization. Errors 
(open boxes) are the combination of statistical (negligible), systematic and normalization uncertainties. 

empirical factor 0.78, at V \v = 900 GeV, n/v = 2.76 TeV and s/s = l TeV [17] are shown. The d/V/dy values 
from PHENIX are reported for particles and antiparticles separately, while the results at LHC energies are 
the average between positively and negatively charged particles, since particle and antiparticle spectra are 
compatible at these energies. We notice that the CMS Collaboration does not include, in the systematic 
uncertainties associated to dN/dy and (pt), the contribution due to the extrapolation to p-\ = 0. For this 
reason, in Fig. 16 and Fig. 17, the ALICE uncertainties are larger than the CMS ones. Similar results 
from the STAR Collaboration [43] are not included, here and in the following plots, since they are 
provided for non-single diffractive events and include contributions of feed-down from weak decays. 

The (K + +K~)/(7r + +7r~) and (p+p)/(7r + +7T~) ratios, as a function of the centre-of-mass energy, are shown 
in the top and bottom panels of Fig. 16, respectively. Results at mid-rapidity from ALICE at sfs = 0.9, 
2.76 [15, 16] and 7 TeV, CMS at ^ = 0.9, 2.76 and 7 TeV [17], PHENIX at ^ = 62.4 and 200 GeV [18] 
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Fig. 16: (K 1 +K )l(n 1 +n~) (top) and (p+p)/(7 T+/r ) (bottom) ratios in pp and pp collisions as a function of the 
collision energy y/s. Errors (open boxes) are the combination of statistical (negligible) and systematic uncertain¬ 
ties. 

and NA49 at y/s = 17.3 GeV [19-21] are displayed. The ratio (p+p)/(7t 1 +7t ) from NA49, calculated 
from the measured particle yields, is not reported because the uncertainty can not be computed from the 
results published in [19-21]. Results in proton-antiproton collisions from E735 at y/s = 0.3, 0.54, 1 and 
1.8 TeV [22, 23] and UA5 at y/s = 0.2, 0.546 and 0.9 TeV [24] are reported, but a direct comparison with 
them is not straightforward due to different baryon number in the initial state. The E735 Collaboration 
provides measurements only for p and not for p yields. Hence the proton-to-pion ratio is computed as 
2p/(7T + +7T - ). In addition, the E735 results for the proton-to-pion ratio are shown in Fig. 16 only for y/s = 
1.8 TeV because at the other energies the p spectra include contributions of feed-down from weak decays 
and arc not directly comparable with the measurements provided by the other experiments. For y/s > 0.9 
TeV, no dependence on the centre-of-mass energy of the (K + +K _ )/(7r + +7r~) and (p+p)/(7T l +K ) ratios 
is observed within uncertainties. 

In Fig. 17 the average transverse momenta (pj) of pions, kaons, and protons, extracted from the sum of 
particle and antiparticle spectra, as a function of the centre-of-mass energy are reported. Results at mid¬ 
rapidity in proton-proton collisions from ALICE at = 0.9, 2.76 [15, 16] and 7 TeV, CMS at /s = 0.9, 
2.76 and 7 TeV [17] and PHENIX at y/s = 62.4 and 200 GeV [18] are shown. In addition measurements 
obtained with E735 at y/s = 0.3, 0.54, 1 and 1.8 TeV [22] and UA5 at y/s - 0.2, 0.546, 0.9 TeV [24] in 
proton-antiproton collisions are also reported. The values of (p/ of p from E735 are not shown since the 
spectra include contributions of feed-down from weak decays and hence are not directly comparable with 
the values provided by the other experiments. A slight increase of (pj) with increasing centre-of-mass 
energy is observed. This rising trend is in particular apparent for /s > 0.9 TeV and it could be related 
to the increasing importance of hard processes at these energies. At /s = 7 TeV, the ALICE and CMS 
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results are consistent within uncertainties except for the proton (pt). This discrepancy is mostly due to 
the difference in the shape of the proton spectra seen in Fig. 14, rather than to the extrapolation to the 
unmeasured p-\ range: a 13% difference is observed on the (pj) values calculated from the ALICE and 
CMS data points in the common p \ range. 



Fig. 17: (Colour online) (pj) as a function of the centre-of-mass energy. Errors (open boxes) are the combination 
of statistical (negligible) and systematic uncertainties. Normalization uncertainties are not included. 


4 Comparison to models 

The comparison between the measured pj spectra of n, K and p and the calculations of QCD-inspired 
Monte-Carlo event generators gives useful information on hadron production mechanisms. Figure 18 
shows the comparison of the measured pion, kaon and proton pj spectra, sum of particles and antipar¬ 
ticles, with two tunes of the PYTHIA6 generator (PYTHIA6-CentralPerugia2011 [26] and PYTHIA6- 
Z2 [27]) 1 , PYTHIA8 tune 4Cx [2, 3], EPOS LHC [4, 5] and PHOJET [6]. 

These event generators are often used and tested to describe hadron collisions at high energies. PYTHIA 
is a general-purpose pQCD-based event generator, which uses a factorized perturbative expansion for 
the hardest parton-parton interaction, combined with parton showers and detailed models of hadroniza- 
tion and multiparton interactions. All presented PYTHIA tunes use a color reconnection mechanism [1] 
which can mimic effects similar to that induced by collective flow in Pb-Pb collisions [45]. In both 
PHOJET and EPOS, which are microscopic models that utilize the color-exchange mechanism of string 
excitation, the hadronic interactions are treated in terms of Reggeon and Pomeron exchanges. 
PYTHIA6-Z2 tune is based on the first measurement of multiplicity distributions in minimum-bias pp 
collisions at sjs = 900 GeV at the LHC. In the CentralPerugia2011 tuning both LEP fragmentation func¬ 
tions and minimun-bias charged particle multiplicity and underlying event data from the LHC are used. 
Both PYTHIA8 and EPOS LHC are tuned to reproduce the existing data available from the LHC (e.g. 
multiplicity and, for EPOS, also identified hadron production up to 1 GeV/c for pions and kaons and up 
to 1.5 GeV/c for protons). The PHOJET parameters are not retuned using the LHC data. 

The measured pion pj spectrum is reproduced by EPOS within 15% over the whole pj range. PYTHIA6- 
Z2, PYTHIA6-CentralPerugia2011 and PYTHIA8 show similar trends. They correctly predict the shapes 
of the pion spectra for pj > 500 MeV/c, overestimating the data by about 10%, 20% and 25%, respec¬ 
tively, while the shapes differ from data for pj < 200 MeV/c (the ratios are not flat) and the yields are 
underestimated by up to 30%. The PHOJET generator does not provide a satisfactory description of the 
measured spectrum shape for any of the particle species. The deviations from the data show a maximum 
for pj ~ 1.2 GeV/c and are more pronounced for kaons and protons than for pions. All the tested Monte- 
Carlo generators underestimate the kaon yield by about 20-30% for pj > 600 MeV/c while for p-\ < 

'The PYTHIA6 tunes are simulated using Rivet [44], a toolkit for validation of Monte-Carlo event generators. 
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Fig. 18: Top panel: Measured pj spectra of pions, kaons and protons, sum of particles and antiparticles, compared 
to PYTHIA6-Z2, PYTHIA6-CentralPerugia2011, PYTHIA8, EPOS LHC and PHOJET Monte-Carlo calculations. 
Statistical (vertical error bars) and systematic (open boxes) uncertainties are reported for the measured spectra. 
Bottom panels: ratios between data and Monte-Carlo calculations. 

400 MeV/c they overestimate the data by up to 30%. A similar deviation is observed by the ALICE col¬ 
laboration also for other strange particle species with a hierarchy depending on the strangeness content 
[46]. The proton yield is well described by EPOS only at low transverse momenta (pj < 1 GeV/c), while 
the generator tends to overestimate the data by up to 30% at higher pj. None of the three PYTHIA tunes 
describes the shape of the proton spectrum in the full pj range. All of them give a reasonable description 
of the yield in the range 1 < pj < 2 GeV/c, but they overestimate the data at lower and higher pj by up 
to 40%. 

The comparison of the pr-dcpcndcnt particle ratios with models allows the hadronization and soft parton 
interaction mechanisms implemented in the event generators to be tested. In the left and right panels of 
Fig. 19, the measured (K 1 +K )/(7i 1 +K~) and (p+p)/(7r + +7r~) ratios as a function of p-\ are compared 
with the same event generators shown in Fig. 18. The measured (K 1 +K ')/(K~+n ) ratio increases from 
0.05 at pj = 0.2 GeV/c up to 0.45 at pj ~ 3 GeV/c with a slope that decreases with increasing pj. All 
the models underestimate the data at high momenta, with EPOS exhibiting the smallest deviation. The 
measured (p+p)/(7T + +7z;~) shows an increase from 0.03 at pj = 0.3 GeV/c up to 0.25 at pj ~ 1.5 GeV/c 
while above this pj it tends to flatten. The data are well described by PYTHIA6-Z2 while PYTHIA6- 
CentralPerugia2011, PHOJET and EPOS show a large deviation at high momenta. PYTHIA8 shows a 
smaller deviation over the whole momentum range even if, as seen in Fig. 18, it overestimates both pion 
and proton spectra. 
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Fig. 19: Measured (K + +K )/(n + +K ) (left) and (p+p)/(7T + +7r ) (right) ratios as a function of pj compared 
to PYTHIA6-Z2, PYTHIA6-CentralPerugia2011, PYTHIA8, EPOS LHC and PHOJET calculations. Statistical 
(vertical error bars) and systematic (open boxes) uncertainties are reported for the measured spectra. 

The comparison between data and Monte-Carlo calculations shows that the tunes of the generators based 
only on few global observables, such as the integrated charged hadron multiplicity, allow only for a 
partial description of the data. The high-precision measurements of identified charged hadron p-\ spectra 
reported here, which cover a wide momentum range in the central rapidity region, give useful information 
for a fine tuning of the Monte-Carlo generators and a better understanding of soft particle production 
mechanisms at LHC energies. 

5 Summary 

A detailed analysis of primary n ± , K*, p and p production in proton-proton collisions at yfs = 7 TeV with 
the ALICE detector has been performed. Particle identification is performed using several techniques 
namely the specific ionization energy loss measured in the ITS and TPC, the time of flight measured 
with the TOF detector, the Cherenkov radiation measured in the HMPID and the kink topology identifi¬ 
cation of the weak decays of charged kaons. The combination of these techniques allows for precision 
measurements of the pj spectra over a wide momentum range: from 0.1 up to 3 GeV/c for pions, from 0.2 
up to 6 GeV/c for kaons and from 0.3 up to 6 GeV/c for protons. A comparison of the ALICE results with 
similar measurements performed by the PHENIX Collaboration at RHIC shows that the /q -intcgratcd 
yields increase with collision energy for all the measured particle species. A slight increase of the (pj) 
with \/s is also observed. This rising trend that becomes apparent at y/s > 0.9 TeV is established by the 
higher y/s LHC data. It could be related to the increasing importance of hard processes at these energies. 
The /q-integrated K/n and pin ratios extend the measurements available at lower collision energies from 
SPS, SppS and RHIC experiments showing a saturation above y/s = 0.9 TeV. Finally, the /q spectra and 
particle ratios have been compared with the calculations of QCD-inspired Monte-Carlo models namely 
PYTHIA6-Z2, PYTHIA6-CentralPerugia2011, PYTHIA8, EPOS LHC and PHOJET. Even though the 
shapes of the spectra are fairly well reproduced by all models (except PHOJET that fails to describe the 
spectrum shape of all the three hadron species), none of them can describe simultaneously the measured 
yields of pions, kaons and protons. These results can be used for a better understanding of the hadron 
production mechanisms in pp interactions at LHC energies and could further constrain the parameters of 
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the models. 
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